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Abstract 

> ' 
O 

■ We propose a systematic method of analyzing pseudopotential transferability 

try 

■ based on linear-response properties of the free atom, including self-consistent 
chemical hardness and polarizability. Our calculation of hardness extends the 
approach of TeterS not only by including self-consistency, but also by general- 
izing to non-diagonal hardness matrices, thereby allowing us to test for trans- 
ferability to non-spherically symmetric environments. We apply the method 
to study the transferability of norm-conserving pseudopotentials for a variety 
of elements in the Periodic Table. We find that the self-consistent corrections 
are frequently significant, and should not be neglected. We prove that the 
partial-core correction improves the pseudopotential hardness of alkali metals 
considerably. We propose a quantity to represent the average hardness error 
and calculate this quantity for many representative elements as a function of 
pseudopotential cutoff radii. We find that the atomic polarizabilities are usu- 
ally well reproduced by the norm-conserving pseudopotentials. Our results 
provide useful guidelines for making optimal choices in the pseudopotential 
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generation procedure. 
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I. INTRODUCTION 



Density-functional calculations performed within the framework of the local-density ap- 
proximation (LDA) have been demonstrated to give accurate predictions of many physical 
properties of solids.! The introduction of the pseudopotential approximation greatly simpli- 
fies electronic-structure calculations by eliminating the need to include atomic core electrons 
and the strong potentials responsible for binding them.! The introduction of norm-conserving 
pseudopotentials by Hamann, Schliiter and Chiang (HSC)@ led to greatly improved control 
of transferability errors, and as a result the pseudopotential approach has since found a wide 
range of applications in molecular and solid-state electronic-structure theory. Nevertheless, 
transferability is still an issue in many calculations, especially when uncomfortably large 
pseudopotential cutoff radii have been dictated by the requirements of a modest plane-wave 
cutoff in the solid-state calculation, and for atoms having shallow core shells. 

A pseudopotential (PSP) is constructed to replace the all-electron (AE) atomic potential 
in such a way that core states are eliminated. The most important measure of a pseudopo- 
tential is its transferability, which characterizes the accuracy with which it mimics the real 
AE atom in different atomic, ionic, molecular, or solid-state environments. Traditionally, the 
transferability of a pseudopotential is characterized by two properties: (i) a comparison of 
the scattering properties of the real and pseudo versions of the free atom or ion, as measured 
by the logarithmic derivative of the wavefunction at some diagnostic radius as a function of 
energy; and (ii) configuration tests, which check if the pseudo eigenvalues and total energies 
track the AE ones for various excited states of the free atom or ion. It is important to note 
that spherical symmetry is implicit in both of these approaches, so that neither is capable 
of giving information about transferability to anisotropic environments. 

Scattering properties are certainly a significant aspect of transferability: poor scattering 
properties are indicative of a poor pseudopotential. A major contribution of HSC was to 
show that the norm-conserving condition automatically implies that not only the logarith- 
mic derivative, but also its energy derivative, is guaranteed to be correct for the PSP at a 
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reference energy, thus insuring that the AE and PSP scattering properties will match closely 
over a large energy range.! Thus norm-conserving potentials tend to have much better trans- 
ferability for most elements. However, the matching of logarithmic derivatives is not always 
sufficient to ensure good transferability. Some potential sources of error which will not show 
up in tests of scattering properties are (i) electrostatic screening effects, and (ii) effects of 
non-linearity of the exchange-correlation energy (important for many alkali- metal elements). 
Errors of the former type are usually easily eliminated by the choice of a "conservative" cut- 
off radius,i@ while those arising from core- valence overlap can largely be corrected by use of 
the frozen-core correction.! Nevertheless, these examples illustrate the dangers of focusing 
on scattering properties alone. 

Configuration tests are certainly useful as a supplementary criterion, but as mentioned 
above, they do not control the quality of the PSP in a nonspherical target environment. 
This will obviously be important for atoms in surface, defective, molecular, or liquid envi- 
ronments, to name just a few. Furthermore, it is difficult to include the configuration tests 
systematically as part of the PSP generation procedure. 

As LDA calculations are pushed in the direction of high accuracy, PSP errors become 
less tolerable, putting tight requirements on transferability. On the other hand, as the calcu- 
lations are pushed to larger system sizes, the increased computational load requires that the 
PSP be as smooth (soft) as possible. This has led to a tremendous effort to optimize PSP 
softness. Unfortunately, transferability and softness are usually contradictory require- 
ments. Especially for first-row elements, attempts to save computational cost frequently 
result in the use of a PSP with an uncomfortably large core cutoff radius. Because of the 
above-mentioned electrostatic screening problems, this strategy may result in a sacrifice 
of transferability which would be difficult to detect using the conventional methods. It is 
therefore of great importance to develop improved measures of transferability which will 
allow for improved control of PSP errors in cases like these. It is especially desirable to 
develop methods which work directly at the atomic level, without the need for painstaking 
comparisons of pseudo and AE results in molecular and solid-state environments. 
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In this paper, we propose to use the linear-response properties of a reference free atom 
or ion as a measure of the transferability of a PSP. We calculate two kinds of linear-response 
properties: a generalized chemical hardness, and the dipole and higher-moment suscepti- 
bilities. Chemical hardness measures the derivatives of electronic eigenvalues with respect 
to changes of occupation. It was recently proposed by TeterB as an important measure of 
PSP transferability, partly based on the idea that the chemical potential and hardness have 
equal roles in determining electron charge transfers. To some extent, the hardness analysis 
is redundant with configuration tests (of which it is a kind of differential version), but it 
is more systematic and can be incorporated into PSP generation more easily.0 The concept 
of hardness is generalized in this paper to include also information about the response to 
non-spherical perturbations. This is in fact very important, and we shall see that the rear- 
rangement of charge in a p-shell may dominate the error in the hardness matrix. We further 
include the self-consistent change of the wavefunctions in our calculations in order to go 
beyond the frozen-wavefunction approximation (FWA) introduced by Teter.0 Finally, the 
dipole (and quadrupole, etc.) susceptibility tests measure the ability of the pseudo atom to 
imitate the correct AE behavior in a local electric field (or field gradient, etc.) which may 
result from an anisotropic solid-state or molecular environment. 

This paper is arranged as follows. In Sec. II we review the concept of chemical hardness 
and introduce the formulation for calculating both spherical and non-spherical hardness. We 
also outline the calculation of the dipole and higher-order polarizabilities. In Sec. Ill, we 
present calculated chemical hardness matrices and polarizabilities for some representative 
atoms chosen from different parts of the Periodic Table. We discuss the general trends in the 
hardness matrix, and the effect of the frozen-core correction. We also proposed a quantity 
based on the hardness matrix to characterize the pseudopotential transferability. Finally, 
we conclude in Sec. IV. 



II. THEORY 
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A. Chemical hardness 

Teter defines the chemical hardness matrix within the LDA as0 

H - 1 m 

" ~ 2 dfidfj ' ' 

where E is the Janak functional^ and fi is the occupation number of the zth state. Since 
the eigenvalue 6j = dE/dfi, we have 

Thus, the hardness matrix measures the first-order change of an energy eigenvalue resulting 
from a first-order variation of an occupation number, while allowing the total number of 
electrons to vary. It consists of two parts: one due to the change of the screening potential 
with variation of the occupation number for fixed wavefunctions, and one arising from re- 
laxation of the wavefunctions. Teter made the approximation of omitting the second term, 
but both are included here. 

Frequently, one is interested only in the case where the occupations of the states compris- 
ing a given shell are kept equal. For example, one may consider an excitation in which one 
transfers an s electron to the p shell, increasing the occupation of each p state by 1/3. This 
insures retention of spherical symmetry of the charge density and potential, and is implicit 
in all of the analysis which is usually carried out with an atomic pseudopotential program. 
In this case, different m components remain degenerate and the treatment is simple .0 In real 
situations (in molecules, at surfaces, etc.), atoms may have very anisotropic environments, 
so that non-spherical changes of electron occupation become important. This prompts us 
to consider also changes of occupation which lead to nonspherical changes of density and 
screening potential. We will use the index L to refer to density or potential changes having 
angular character Ylm(^)- Thus, we shall not restrict ourselves to spherically symmetric 
(L = 0) perturbations, but will consider the general L ^ case. 

For this purpose, it is useful to generalize from the concept of an occupation number 
fi to the concept of an "occupation matrix" or "density matrix" fy. This generalization, 
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previously introduced in other contexts,E3'L3 makes the hardness analysis more complete. It 
is described in Subsec. IIC below. While fy is diagonal in the atomic ground state, or in a 
basis of energy eigenstates of a perturbed system, it may be non-diagonal in a more general 
representation of a perturbed system such as an atom in a defective environment. 

In our calculations, we use first-order density-functional perturbation theory in the frame- 
work of LDA, following the scheme formulated by Mahan and SubbaswanryB In the remain- 
der of this section, we first give a detailed formulation of the calculation of the hardness 
matrix elements associated with conventional diagonal occupation number changes, and 
discuss the extension beyond the frozen- wavefunction approximation (FWA). Next, the gen- 
eralization to non-diagonal occupation changes will be presented. Finally, we sketch the 
calculation of the dipole and higher susceptibilities, which is straightforward after the ma- 
chinery needed to calculate the hardness elements has been set up. 



B. Hardness for diagonal occupation 

Within the framework of LDA, application of the Hellmann-Feynman theorem to Eq. 
(|) yields 

l OJj 

= , (3) 

2 mi dp dfj m/ ' 1 ; 

where Vh xc is the Hartree and exchange-correlation potential. With atomic quantum num- 
bers, this becomes 

Hnlm =\ I I dr dr' n nlm (r) w hxc (r,r') - . (4) 

n'l'm' I J J Otn'Vm' 



Here, 



n nlm {r) = R 2 nl (r) |l^ m (fi)| 2 , (5) 

, f\ dVhxc(r) 
Whxc (r,r)= , (6) 

occ 

p(r) = fmm n n im{r) ■ (7) 

nlm 
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The density variation due to the variation of occupation numbers consists of two terms: 

j-, — = n n i m (r) + An nlm (r) . (8) 

The first term arises from the explicit dependence of density on the occupation numbers, 
while the second involves relaxation of the wavefunction with changes of occupation. We 
will refer to the neglect of An n i m as the "frozen wavefunction approximation" (FWA), while 
the effect of the An n ; m term will be referred to as the "self-consistency" (SC) correction. 

The FWA part of the hardness matrix is relatively easy to calculate, and it has previ- 
ously been done for almost all atoms in the Periodic Tabled The SC correction is treated 
using density-functional linear response theory, regarding the change of f n i m as an external 
perturbation. 

We start with the calculation of the FWA hardness, which may be written 



rrFWA 
tl. 



n'l'm 

where 



nlm 



- / drV n i 

in 

(r) n n > v 

ml V / > (9) 



V n i m {r) = J dr' w hxc (r, r') n nlm (r') . (10) 
The kernel tfhxc can be decomposed as 

w hxc (r, r') = ^ «ffi(r, r') Y£ M (Q) Y LM (Q') , (11) 

LM 

where 

(L) ( f\ 8ti r <L SV XC 5{r-r') 

and r < and r > are the smaller and larger of r and r' respectively. Note that the exchange- 
correlation term, being local, is independent of L. The functions n n i m and V n i m can also be 
expressed in spherical harmonics as 

21 

n nlm {r) = C^{l,m) n nl {r) Y L)0 {tt) , (13) 

L=0 



21 

W) = E ^ (L) (^,H v$\r) F Li0 (O) , (14) 

L=0 

where the are Gaunt coefficients, reflecting the presence of squared spherical harmonics 
in n n i m (r). Here and in the remainder of this subsection, the sums are over even L only. 
Defining 
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J drdr' r 2 r' 2 n n i(r)w\^ c (r,r')n nlu (r') , (15) 



the FWA hardness can be expressed as 



1 2/ m i n 

^mm = ? E C (L) (/, m) CW(l'm') ^rm/u'L , (16) 



where l m i n is the smaller of I and Z'. 

Of course, the presence of V n i m (r) induces a change An„[ m (r) in the charge density. The 
SC correction to the hardness can be calculated by including this change self-consistently, 
treating V n \ m as a bare potential perturbation. We can write 

A» nim (r) = J dr' X (r, r') V nlm {v') (17) 

where \ is the linear susceptibility (to be discussed below). According to a theorem by 
Eaves and Epstein,0 the induced charge density for a closed-shell atom has the same angular 
character as the perturbing potential. Within LDA (with spin-polarization neglected) the 
ground state is always isotropic (equal population of each m character of any given shell), 
so this theorem applies. Thus, to linear order we can write 

Ang^r) = f dr' r' 2 x {L \r,r') v%\r') (18) 



where ( r > r ') is the linear susceptibility in radial coordinates for perturbations of angular 
character L, and 

21 

An nim (r) = CM(l,m) Anjf } (r) Y Lt0 (n) , (19) 

L=0 



in analogy with Eq. flTq). After defining 

Uffiwr. = \ J J dr dr' r 2 r' 2 v£> (r) x (L) (r, r') (0 , (20) 
the SC correction to the hardness matrix becomes 

AH n im =- f [ drdr'V n i m (r)x{r,r')V n nr m '{r') 

n'l'm' 2 J J 

= \ 2l f:C (L \lm)C^\l'm')Ahl% WL . (21) 

Z L=0 

The determination of the linear susceptibility x(r, r') follows closely the modified Stern- 
heimer approach discussed by Mahan and Subbaswamy.0 In what follows, we consider the 
response An^(r) = An^(r)Yko(fi) to a general perturbation v^ l \t) = v^ L ' (r)Yi (Q) of 
angular character L. (We have v = v n i and An = An n i in the immediate context, where 
nl are indices of the state whose occupation is being varied. However, the discussion given 
below is general, and the indices nlm will henceforth be taken to refer to arbitrary occupied 
wavefunctions which respond to the perturbation.) We do not actually need to calculate \ 
itself; instead, it is sufficient to specify an iterative algorithm for calculating its action 

Afi (i) (r) = / dr' r' 2 X {L) (r, r') v (L) (r) (22) 



upon the arbitrary perturbation. 

The procedure is as follows. For the moment, assume that the first-order density change 
An^(r) is known; then the corresponding change Av^(r) in the screened Kohn-Sham 
potential is given by 

Aug 5 (r) = v {L) (r) + J dr' r' 2 w$ c (r, r') An {L) (/) . (23) 

This induces a first-order change in each Kohn-Sham wavefunction satisfying 

(Ho - e nl ) Aip nlm = - (A^f' - Ae nlm ) tp nlm , (24) 



where Ho,e n i, and ^ n j m (r) are the unperturbed Kohn-Sham Hamiltonian, eigenvalue, and 
eigenfunction, respectively, and Ae n i m and Ai/; n i m (r) are the corresponding first-order 
changes. Decomposing Eq. (|24]) into angular and radial parts, one finds 
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l+L 

A^w(r)= £ ^(^0;/m,/ , m)A J R n ^(r)^ m (fi) (25) 

Z'=]Z-£| 

where the C(LO;lm,l'm) are Clebsh-Gordon coefficients and the radial functions AR n iv{r) 
are the solutions of the radial part of Eq. ( p|) . 



■^-2 + 2 + M on (r) + V hxc (r) - e nl 



AR nlv {r) = [Ae nlrn - Av sci (r)] R nl (r) . (26) 
This inhomogeneous equation is solved numerically on a radial mesh following Ref. in. 



Finally, the change of density resulting from Eq. fl2"5]) is 



An^(r) = £ ^TTT) AjRS(r) ^(WW) > ( 27 ) 



where f n i = J2 m fnim, and the geometric coefficients D(W, L) are given on p. 55 of Ref. [13. 
Iterative solution of Eqs. fl23|), (|26|) , and (^) thus gives the self-consistently screened density 
change An^'(r) resulting from the perturbation v^ L \r) in Eq. (E3). 



C. Non-diagonal hardness 

If one wants to express any atomic charge density in terms of (both filled and empty) 
Kohn-Sham orbitals of a reference ground-state atom, then in general the occupation num- 
bers turn out to be non-diagonal, so that the total charge density should be expressed as 

"00 =E &V£(r)V>i(*) (28) 

where = /, <5 i:; - only in the ground state. The most general result of a perturbation is thus 
obtained by allowing for non-diagonal terms to exist JllEl and by taking this into account 
we arrive at a more general form of the hardness matrix. This provides a natural way to 
test the transferability of a PSP to nonspherical environments. 

Instead of working in the explicit atomic representation /y = f n im,n'i'm'j we find it 
convenient to work in a representation (nn'll'LM) in which LM are labels of total angular 
momentum, and L = \l — \ l — l'\ + 2, I + V following the usual angular-momentum 
addition rules. Introducing also the condensed notation a = nn'll', we thus have 
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n(r)= J2 faLMn a (r)Y LM (n) (29) 

aLM 

where 

faLM = C(LM; Im, l'vn!)f n[m , n , Vm , (30) 

ram! 

with C(LM;lm,l'm') again being the Clebsh-Gordon coefficients, and n a (r) = 
R n i{r)R n iii{r). In this notation, it is natural to introduce the generalized Kohn-Sham eigen- 
values 

dE 

taLM — TT? W-U 

OjaLM 

and the generalized hardness matrix 

tlaf3LL>MM' — W7 ■ 

* OjaLM OjpL'M' 

Since our reference unperturbed configuration is spherically symmetric, e a LM is only 
nonzero for L = M = and n — n', I — I'. Spherical symmetry also implies that the 
hardness matrix reduces to 

Ha(3LL'MM> = H a f3L &LL' &M-M' ■ (33) 

The FWA hardness becomes just 



H F J L A = tj / jdrdr' r 2 r' 2 n Q (r) «fi(r, r') n,(r') , (34) 



1 

2 

where the quantity w^ c is as defined in Eq. ([12]) (now generalized to odd as well as even L) . 
The calculation of the SC correction proceeds along similar lines as for the diagonal case. 
The perturbing potential is generalized to 

V aLM (r) = J dr' r' 2 w^(r, r') n a (r') Y LM {9) 

= v a L) (r)Y LM (Sl) (35) 

and the self-consistent correction is given by 
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AH %3L = \ J Jdrdr' r 2 r' 2 (r) X {L) {r, r') vf ] (r') . (36) 

In the special case a = nnll and j3 = n'n'l'l', the quantities H^m^ and H^9 L reduce 
to h™wL and hf^i H i L , respectively, as defined in the previous subsection. Thus, all of the 
diagonal hardness matrix elements are contained as special cases of the generalized non- 
diagonal ones introduced here. Moreover, the original diagonal formulation only covers 
variations of the screening potential or density of even L (monopole, quadrupole, etc.), 
whereas the generalized formulation is capable of treating variations of any angular character. 
In fact, perturbations of dipole (L = 1) character are likely to be the most important 
nonspherical perturbations in many molecular and solid-state environments, especially at 
surfaces and other defects where inversion symmetry is lacking. Therefore, in what follows 
we will concentrate on comparisons of the AE and PSP non-diagonal hardness elements 
H a /3L, with special emphasis on the L = and L = 1 cases. 



D. Polar izability 

The polarizability of an atom measures its response to an external electric field. The 
dipole (L — 1), quadrupole (L — 2), and higher (L > 2) polarizabilities are defined as the 
derivatives of the L'th induced charge moment with respect to an electrostatic potential 
of form r L Yio(Q). For good transferability, it is important that the pseudo-atom have 
polarizabilities similar to those of the all-electron atom. We can expect the lower-moment 
polarizabilities to be more important, so we focus on the dipole and quadrupole susceptibility 
in what follows. Tests of the non-selfconsistent polarizability of HSC pseudopotentials have 
previously been performed for a large number of closed-shells atoms and ions by Bachelet 
et a/.0 Here, we extend the tests to other atoms and also include the SC correction. 

The formulation and calculation of the polarizability is straightforwardBEl using the 
machinery developed in the previous subsections. The perturbing potential is taken to have 
the form 
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V L (r) = r L Y L0 {Q) 



(37) 



and the linearly induced density change is 



An {L \r) = 



dr' r' L+2 X (L) (r,r') . 



(38) 



This is evaluated using the same iterative procedure given previously in Eqs. (p3|427p . The 
polarizability in angular channel L is then 



In this part, we present our calculated hardness matrix elements for a set of representative 
atoms. We begin with all-electron atoms, surveying the characteristics of the hardness matrix 
and the basic trends as a function of position in the Periodic Table. We concentrate first on 
argon, making the comparison between the AE and PSP hardness matrix, discussing sources 
of error, and introducing our format for presenting results in a systematic fashion. Next, 
we discuss trends as one goes across the Periodic Table, and evaluate the importance of the 
Louie- Froyen-Cohen (LFC)I semi-core correction and the effects of varying the core radius. 
We will then propose a method to extract the most important information from the large 
number of hardness elements. Finally, the calculated AE and PSP polarizabilities will be 
presented at the end of this section. 

The results reported here are restricted to potentials of the HSC type. We believe 
that these results can be taken as indicative for the whole class of norm-conserving HSC- 
like PSPs,@@0'iiHii provided comparable cutoff radii are chosen. For Kleinman-Bylander,0 
ultrasoft,0 or other approaches to PSP construction which deviate significantly from the 
original HSC method, addition investigation may be appropriate. 




(39) 



with pi and pi being the dipole and quadrupole susceptibility, respectively. 



III. RESULTS 
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A. Hardness matrix for all-electron atoms 



Before studying the transferability of PSPs with reference to the hardness matrix ele- 
ments Hij, we first need to understand the behavior of for all-electron atoms in order to 
obtain some physical intuition. In Table |, we show some matrix elements of the generalized 
hardness for three characteristic AE atoms in the Periodic Table. Silicon and argon are taken 
in the ground-state neutral configuration, while for sodium we choose an ionized valence con- 
figuration 3s a5 3p a25 in order to ensure that the p electron is bound. The total self-consistent 
hardness H total = H FWA + H sc is broken down into frozen- wavefunction-approximation and 
self-consistent-correction pieces as discussed in the previous section. For some elements, the 
FWA contribution is further separated into Hartree (h) and exchange-correlation (xc) con- 
tributions, H FWA = H h + H xc . The matrix elements can be identified by quantum numbers 
a = riin^lil'i, (3 = n^n!^^-, anc ^ ^' Since the principle quantum numbers n are obvious in 
most cases, we will usually omit them and write the indices simply as I2I2', L. 

From Table I, we see that H FWA typically dominates, and the self-consistent correction 
H sc is only about 20 — 30% of H FWA . Nevertheless, this is clearly large enough that complete 
neglect of H sc is unjustified. It is also obvious from the table that H FWA are mostly positive, 
while the H sc are negative. For the diagonal elements (e.g., excluding {sp, sp; L = 0}), this 
can be understood as follows. From Eq. fl3~3|), the Hartree part of H FWA can be seen to have 
the form of the Coulomb self-energy of a particular charge distribution n a (r)y^j\/(fi), which 
must be positive. On the other hand, the exchange-correlation contribution is negative 
because SV xc /5n in Eq. (fl2]) is negative. Typically, the Hartree term dominates and H FWA 
is positive. Regarding H sc , note that Eq. (|36[) can be interpreted as evaluating the second- 
order energy change of the system when the external potential of Eq. (^) is applied; the 
total energy must go down when the wavefunctions relax, so that H sc must be negative. 

Going from Na to Ar, we find that H total increases strongly, primarily because of the 
Hartree contribution H h to H FWA . This is in agreement with our intuition, since the wave- 
functions become much more localized as one moves from left to right across the Periodic 
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Table, causing H h to increase sharply. In this sense, we can refer to atoms (like Ar) on the 
right side of the Periodic Table as "strongly electrostatic atoms" or "hard atoms," while 
those on the left side (like K) can be termed "soft." 

Our results also indicate that the matrix elements of H total for L = are significantly 
larger than those for L > 0. This implies that the diagonal hardness elements defined in Eqs. 



(|i6|) and ([H]) are dominated by the spherically symmetric part of the response. A closer 
inspection reveals that H h decreases strongly with increasing L, while H xc is smaller and, 
being a local operator, independent of L. This gives rise to an overall small H total for high 
L. In a few cases where the valence wavefunction was very weakly bound and delocalized, 
we have found that H xc can even be larger in magnitude than H h , resulting in a small but 
negative H total . We regard this as an unphysical result which reflects the overestimate of 
exchange-correlation effects by the LDA in the low-density tail of the atom. 



B. Argon: AE and PSP hardness 

Having gained some understanding of the AE hardness matrix, we turn now to a com- 
parison of the PSP hardness elements with the corresponding AE ones. We begin with 
argon. In Table || we list values for some important hardness matrix elements calculated for 
the AE and PSP Ar + ion in configuration s l 2 p 5 ' r dP' 1 . The HSC PSP was generated in this 
configuration, using a core radius r c ~ 1 a.u. We also show the relative errors for another 
valence configuration, s 1 p 5 ' 5 d ' 5 . 

Comparing the AE and PSP results, we find very good agreement for the L = matrix 
elements. Norm conservation imposes the constraint that the L = component of the 
electrostatic potential in the PSP case should match the AE one outside the core region. 
Consequently, for small core radii the differences between AE and PSP values of H h are 
essentially confined to L > moments. Actually, we find that the errors in H h are relatively 
insensitive to a modest increase of core radius for hard atoms like Ar. Because of this 
constraint, while the Hartree contribution is large in magnitude, it may only incur a small 
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error in hardness. For Ar, the H xc are quite small and their contribution to the error is 
not very significant. So, for Ar, the agreement between AE and PSP hardness elements 
is excellent in the spherically-symmetric (L = 0) channel. The relative errors for L > 
channels are larger (~ 20% for L = 1 and ~ 4% for L — 2). However, since their absolute 
magnitudes are small, they are not as important. As a result, the overall agreement between 
AE and PSP hardness is very good. We also checked that changing the testing configuration 
affects the results only very slightly. 

In order to avoid presenting numerous cumbersome tables, we have converted the infor- 
mation into the form of a bar chart as shown in Fig. p]. The heights of the bars represent 
the values of the corresponding hardness matrix elements. The columns indicate different 
contributions to a given hardness matrix element, whose indices are labeled by the row. The 
results for all-electron and pseudo atoms are placed side-by-side to facilitate comparison; 
hollow bars indicate AE results while dashed bars represent PSP results. The values for 
L > elements, being small, are magnified in the diagram. We believe the discussion is eas- 
ier to follow by viewing such diagrams, so all subsequent hardness results will be presented 
in this way. 

C. From Ar to K 

Argon is a rare-gas element. To explore the general trends of the hardness matrix along 
the Periodic Table, we further calculated the hardness matrix for atoms having a wide range 
of properties. 

We start with Si. To study the effect of core radius r c on the quality of PSP generated, 
we plot in Fig. the results for both small and large r c indicated by dashed and shaded bars 
respectively (values for r c are indicated in the captions). The AE results are still plotted 
with hollow bars. These results are for Si + in configuration s 2 p 1 d°. For small r c , the PSP 
hardness elements for Si agree with the AE results very well (the biggest errors occur for 
L = 1 elements as for Ar). With the exception of {dd,dd;L = 0} elements, a worsening 
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of the agreement in total hardness is evident in the L = channel for larger r c . (The d 
wavefunction is nodeless and very delocalized, and thus insensitive to changes in core radius.) 

One expects that increasing r c should always make the PSP less transferable, but the 
sensitivity can be different for different elements. In Fig. 0, we show a similar diagram for 
oxygen, again using hashed and shaded bars to represent hardness elements for a PSP with 
small and large r c . We find that the effect of increasing r c is more dramatic for 0. A more 
complete picture of the effect of r c will be presented later. 

We next focus on some cases to characterize the role of the LFC correction.! In the 
following figures, we use the shaded bars to represent PSP results with such LFC correction. 
The results we presented are for K +0,25 in s °- 25 p - 25 ^ - 25 (Fig. Ti +0 75 in configuration 
s y.25 d 2 ( Fig> |) ; c +o.4 in s i.y.4 d o.3 ( Fig _ |^ and Ga +o.5 in s 2 p o. 5rf o ( Fig _ ^ Starting 

with K we find as expected that the LFC greatly reduces the error due to the H xc contri- 
bution to the hardness matrix elements even for the L = channel. The LFC-corrected 
hardness matrices are in good agreement with the AE results. For Ti, the LFC successfully 
corrected the noticeable mismatch between AE and no-LFC results for H xc contributions 
to the {dd, dd; L = 0} and {dd, dd; L = 2} elements. (The no-LFC total hardness in the 
{dd, dd; L = 0} happens to match the AE one rather closely, but only because of a fortuitous 
cancellation of errors.) For C and Ga, the effect is small except for some high-L channels 
(i.e. {pp,pp; L = 2} elements), and we regard the LFC as less necessary for these elements. 

The hardness results obviously depend on details of the PSP construction, but we note 
the following general trends. In going from left to right across a row of the Periodic Table, 
an increasing atomic number tends to localize the core density closer to the nucleus and 
to reduce the size of the core. Consequently, the overlap between core and valence charge 
densities, which is the source of the nonlinearity in the exchange- correlation potential, gets 
smaller and smaller. Such an overlap, when significant, is largely responsible for the errors 
in the self-consistent contribution to hardness as well, i.e., for a poor description of the rear- 
rangement of the pseudo wavefunctions. Provided that one takes small core radii to minimize 
errors in the Hartree contributions, the degree of core- valence overlap almost entirely deter- 
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mines the PSP transferability. Rare gas atoms thus have maximum transferability without 
LFC, while alkali atoms need the LFC correction the most. 

Moving along the columns of the Periodic Table, the trends are naturally much weaker. 
The net positive charge seen by the valence electrons doesn't change, and the localization 
of the core and valence charge densities changes only marginally. However, a slight increase 
of overlap occurs, along with a corresponding loss of transferability of non-LFC PSP, as one 
goes down the columns. 

From our hardness matrix results, we can investigate whether there is a systematic way 
to improve the PSP transferability. The Hartree contribution can be improved by imposing 
additional conditions on the pseudo wavefunctions, e.g., matching of the valence electrostatic 
potentials for higher-order multiple moments.il However, it is not clear whether the gain in 
transferability would justify the drawbacks of imposing additional constraints. As regards 
the exchange-correlation terms, while some other approaches have been tried,0 the LFC 
correction seems to be the simplest and most efficient method. 



We have shown that it is useful to characterize the transferability of a PSP in terms of its 
hardness matrix. However, there is so much information contained in the numerous matrix 
elements of the hardness matrix that it becomes difficult to decide whether a particular 
PSP shows "good" or "poor" transferability. Thus, it is desirable to define a single quantity 
that can be used to represent approximately the overall transferability of the PSP. There 
is certainly no unique way to do this, since the importance of different matrix elements 
depends on the target application. Nevertheless, we propose one such definition, which at 
least can be used as a starting point. 

We define an average hardness error X as follows: 



D. Average hardness errors 




(40) 



a/3L 
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Here AH is the difference between total AE and PSP hardness matrix elements, and w Q /5L 
is a weight to be defined shortly. Thus, X is just a weighted RMS average of the errors in 
the hardness matrix elements. 

To fix the weights w a /3L, we have adopted the following philosophy. We want X to repre- 
sent an average total-energy error which would occur as the PSP atom is transferred to an 
ensemble of target environments, where the distribution of target environments is charac- 
terized by specifying the average occupation Ni and the typical fluctuation in occupation rji, 
for each electron shell. In the spirit of the hardness approach, we can estimate the change 
in total energy of the atom as it is inserted into a given environment as 

AE = H *(3L Sf aL M SfflLM • (41) 

a/3LM 

If each contribution were statistically independent, one would have 

( AEf ~ (2L + 1) E HIpl SfaLM 6fi LM ■ (42) 

a/3L 

Making the additional rough approximation that Sf 2 LM oc r]i a r]i' a and 6fp LM oc r]i r]i' (where 
a = n a n' a l a l' a and the n subscripts are suppressed) and replacing (H a pi) 2 by the AE vs. 
PSP error (AH a p L ) 2 , we arrive at the right-hand side of Eq. (|40|) with 

w apL = (2L + 1) 7]i a T] Va Tjip Vi' , (43) 

where it remains to fix the occupation fluctuation rji. To make things simple, we assume 
that the fluctuation rji is a function only of the average occupation iVj and the maximum 
occupation 2(2/ + 1) of the shell. We choose the form 

Vl = (21 + 1) phtl-fi) (44) 

where f] is the fractional occupancy, fi = Ni/2(2l + 1). The first term makes the fluctuation 
proportional to the number of electrons which could be accommodated in the shell, and the 
second forces the fluctuation to zero for either a completely filled or a completely empty 
shell in a manner which respects electron-hole symmetry. 
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Thus, we have defined X through Eqs. (fRJ), (|43|) and fl44]), in such a way that it depends 
only on a specification of the reference electronic configuration (the Ni values). Our definition 
puts no weight on completely filled or completely empty shells, and heavily weights partially- 
filled shells. 

We would be the first to admit that the choices above are largely arbitrary, but we believe 
they are reasonable ones, and we proceed to use this measure to study the effect of variations 
in core radius upon PSP transferability. In Figs. |8| and || we show the calculated average 
hardness error X for a set of six elements as a function of r c = (r cs + r cp )/2. Here, r cp — r cs 
is kept constant. We do not change r c ^, since it does not affect X appreciably for the atoms 
studied here (with no d electrons inside the core). The hardness error X (and the hardness 
itself) is much greater for first-row elements C and O, so they are plotted on a different 
scale. It can be seen that the behavior of X differs considerably between elements. For K, 
X is very insensitive to r c , and there is thus wide flexibility in the choice of an appropriate 
r c . For O, C, and Ga, X increases in an approximately linear fashion as r c is increased from 
1.5 to 3.5 au, while X increases more rapidly for Si and Ar as r c is increased. 



E. Polarizability 



In Table |HI[ we show our calculated dipole (L = 1) and quadrupole (L = 2) polarizability 
for some AE and PSP atoms. Our AE results for rare-gas and closed-shell ions agree very 
well with previous calculations.0i3Il The results for the dipole polarizability are also in 
good agreement with experiment. (Experimental values for higher- moment polarizability do 
not appear to be available.) For example, for the K + ion, we find the dipole susceptibility to 
be 5.74 a.u., compared to 5.86 a.u. and 5.47 a.u. from previous theoryEl and experiment^, 
respectively. Note that most of the results reported below are for open-shell atoms or ions. 
It should be emphasized that in these cases, our results are a theoretical fiction in that 
we assume symmetrized occupations (e.g., s 2 p 2 J 3 p 2 J 3 p 2 J 3 for C) which have little relation 
to the real atomic ground state. Nevertheless, we believe it is meaningful to compare AE 
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and PSP polarizability calculated in this way as a means of testing the transferability of a 
PSP. We tend to prefer tests on ionized configurations (e.g., Si + instead of Si) because we 
have found that shallow orbitals in neutral open-shell atoms sometimes give such enormous 
contributions to the polarizability that comparison becomes difficult. All PSPs are built 
choosing small core radii (e.g., r cs = 1.1 a.u. for Si + ), and for K the LFC correction was 
used. 

We first consider the all-electron results. To test the effect of self-consistent screening, 
we report both frozen-wavefunction approximation (FWA) and self-consistent (SC) polar- 
izabilities. All pseudopotentials are built choosing small core radii, and for potassium, the 
LFC correction is used. In all cases, the core polarizability is included in the all-electron 
value. As a example, we show the core contribution of K + to its all-electron polarizability. 

Looking at Table |TJ, it is evident that only the dipole polarizability is strongly affected 
by self-consistent screening. The screening reduces the dipole polarizability by around 40%, 
while the quadrupole susceptibility is typically reduced by only about 3%. Exceptions to 
this pattern occur for some highly polarizable atoms like potassium, for which the screening 
correction is still sizeable in the quadrupole channel. 

Three factors contribute to the difference between PSP and AE polarizability: (i) the 
core contribution; (ii) the difference between unperturbed pseudo and AE wavefunctions 
inside the core region; and (iii) differences in the first-order changes in the valence wave- 
functions. Regarding (iii), the wavefunction changes are determined in part by admixture 
of angular-momentum components higher than those which are present in the unperturbed 
reference configuration. For these components (typically I > 3), no norm-conservation or 
tail-matching conditions were imposed. Because the PSP usually contains no nonlocal pro- 
jectors for large I, these wavefunctions feel only the local potential, which is usually set in 
a very arbitrary manner. This appears to be the most significant source of error in PSP 
polarizability. For example, we have calculated the PSP polarizability of Ar both with and 
without an / component in the nonlocal projector. We find that the calculated values for 
both the dipole and quadrupole polarizabilities are about 15% too small when the / com- 
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ponent of the projector is omitted. This effect was already noted in Ref. [17|, where many 
other examples can be found. All the results in Table [TTI| are obtained using a PSP with a 
complete projector (up to / = 3). 

The results in the Table indicate that while the screening correction generally improves 
the agreement between AE and PSP for the dipole susceptibility, the error in the quadrupole 
susceptibility is almost unaffected. Generally, we find a very good agreement between the 
self-consistent AE and PSP results. However, it should be noted that all results shown in the 
Table were obtained using a PSP generated from the same configuration for which the polar- 
izability calculation was made. If we change the PSP reference configuration considerably, 
larger changes in the calculated pseudo-polarizability may occur. 



IV. CONCLUSION 

We present a systematic method for characterizing the transferability of pseudopoten- 
tials using their linear-response properties, specifically their generalized chemical hardness 
and dipole and quadrupole polarizabilities. The hardness measures the ability of the PSP 
to resemble the all-electron atom in different atomic environments, including non-spherical 
ones, while the polarizability reflect the response of the PSP atom to external fields. When 
used together with conventional criteria such as norm-conservation and matching of eigen- 
values and logarithmic derivatives, this approach allows a rather complete characterization 
of PSP transferability. 

We have applied the method to study the behavior of Hamann-Schluter-Chiang pseu- 
dopotentials for many atoms in the Periodic Table. As expected, the calculated hardness 
matrix indicates that the transferability deteriorates as the core radius is increased. For 
some elements with relatively delocalized cores, we find strong evidence for the importance 
of including the Louie-Froyen-Cohen semi-core correction. We propose a method for reduc- 
ing the large amount of information contained in the hardness matrix to a single number. 
We suggest that this quantity be monitored or included in the fitting procedure when gener- 
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ating pseudopotentials, in order to achieve the desired properties (e.g., optimal smoothness) 
without sacrificing transferability. 
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TABLES 

TABLE I. Calculated all-electron hardness matrix elements of three representative atoms Na, 
Si, and Ar. Both frozen-wavefunction approximation (FWA) and self-consistent (SC) contributions 
to the total hardness elements are listed. For some matrix elements, the FWA hardness is broken 
down into Hartree (h) and exchange-correlation (xc) contributions. 

Na(3s a5 3p - 25 ) Si(3s 2 3p 2 ) Ar(3s 2 3p 6 ) 

ss, ss; L = 

FWA 0.1906 0.3923 0.6142 



SC 


-0.0166 


-0.0931 


-0.1684 


total 


0.1741 


0.2992 


0.4457 


ss,pp; L = 








FWA 


0.1562 


0.3437 


0.5668 


SC 


-0.0144 


-0.0737 


-0.1444 


total 


0.1417 


0.2700 


0.4224 


pp,pp; L = 








FWA(h) 


0.1625 


0.3260 


0.5465 


FWA(xc) 


-0.0374 


-0.0197 


-0.0197 


FWA 


0.1251 


0.3063 


0.5268 


SC 


-0.0131 


-0.0586 


-0.1240 


total 


0.1120 


0.2477 


0.4028 


sp, sp;L = l 








FWA(h) 


0.0367 


0.0766 


0.1275 


FWA(xc) 


-0.0319 


-0.0203 


-0.0212 


FWA 


0.0048 


0.0563 


0.1063 


SC 


-0.0014 


-0.0306 


-0.0560 


total 


0.0033 


0.0257 


0.0504 



pp, pp; L = 2 
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FWA(h) 

FWA(xc) 

FWA 

SC 

total 



0.0164 
-0.0374 
-0.0210 
-0.0031 
-0.0242 



0.0321 
-0.0197 

0.0124 
-0.0018 

0.0106 



0.0536 
-0.0197 

0.0339 
-0.0068 

0.0272 
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TABLE II. All-electron (AE) and HSC pseudo-potential (PSP) hardness matrix elements for 
Ar. Two different electronic configurations are considered. The error is the percentage difference 
between PSP and AE hardness. Total hardness (total) is decomposed into frozen-wavefunction 
approximation (FWA) and self-consistent (SC) contributions, while the FWA is further decomposed 
into Hartree (h) and exchange-correlation (xc) contributions. 







1 2 5 7 j0 1 

s p d 




1 5 5 j0 5 

s p d 




i\Hj i rvy ) 




error ( /o J 


ei i ui ^ /o J 


ss, ss', L = 










b WA(h) 


0.6588 


0.6562 


n on 

0.39 


0.41 


b WA^xcj 


-0.02/3 


n none 
-0.0308 


12.98 


1 9 no 
13. U9 


b WA 


n £?91 c 
0.6315 


0.6254 


n n^7 

0.97 


1.01 


SC 


-0.1421 


-0.1367 


3.80 


3.26 


total 


n /i qo/1 
U.4oy4 


U.4oo f 


n i r 
U.lO 


U.uY 


pp,pp; L = 










r WA^nJ 


0.58 I 6 


n c c /I c 
0.5845 


n c q 
0.53 


n c /i 
0.54 


b WA(xcJ 


n noon 
-0.023U 


n no ix /i 
-0.0254 


1 n o /i 
10.24 


in ci 
1U.51 


t WA 


0.5646 


n c cm 

0.5591 


n n^7 

0.97 


1 nn 
1.0U 


SC 


-0.1139 


-0.1092 


4.17 


3.44 


rorai 


n a <^c\7 

U.40U / 


n a ^nn 

U.40UU 


U. ID 


nil 

U. 1 1 


dd, dd; L = 










FWA(h) 


0.3237 


0.3234 


0.10 


0.12 


FWA(xc) 


-0.0341 


-0.0344 


0.81 


1.81 


FWA 


0.2897 


0.2890 


0.24 


0.29 


SC 


-0.0380 


-0.0375 


1.24 


1.37 


total 


0.2517 


0.2515 


0.09 


0.05 


sp, sp; L = 1 










FWA(h) 


0.1365 


0.1417 


3.84 


3.91 
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FWA(xc) 


-0.0244 


-0.0274 


12.12 


12.25 


FWA 


0.1120 


0.1144 


2.13 


2.05 


sc 


-0.0475 


-0.0342 


27.91 


27.45 


total 


0.0645 


0.0801 


24.25 


21.17 


sp, pd; L = 1 










FWA(h) 


0.0935 


0.0987 


5.56 


5.87 


FWA(xc) 


-0.0160 


-0.0178 


11.12 


12.03 


FWA 


0.0775 


0.0809 


4.41 


4.58 


SC 


-0.0331 


-0.0254 


23.23 


22.75 


total 


0.0443 


0.0555 


25.30 


22.54 


pp,pp; L = 2 










FWA(h) 


0.0586 


0.0606 


3.33 


3.31 


FWA(xc) 


-0.0230 


-0.0254 


10.24 


10.51 


FWA 


0.0356 


0.0352 


1.13 


1.43 


SC 


-0.0055 


-0.0037 


32.48 


31.92 


total 


0.0301 


0.0315 


4.60 


3.87 


dd, dd; L = 2 










FWA(h) 


0.0283 


0.0286 


1.03 


1.09 


FWA(xc) 


-0.0341 


-0.0344 


0.81 


1.81 


FWA 


-0.0057 


-0.0058 


1.49 


0.96 


SC 


-0.0026 


-0.0022 


15.29 


23.91 


total 


-0.0083 


-0.0080 


3.76 


8.88 
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TABLE III. Comparison between all-electron (AE) and pseudopotential (PSP) dipole and 
quadrupole polarizability, in atomic units (e 2 = 2), for selected ions. 



Frozen-Wave Approximation Self-Consistent Results 





AE 


PSP 


error (%) 


AE 


PSP 


error (%) 


Dipole 














K+ 


8.89 


— 


— 


5.74 


— 


— 


K +0.3 (s 0. 7) 


181.94 


175.76 


3.4 


165.89 


164.83 


0.6 


C+ (sV) 


10.04 


10.16 


1.2 


6.04 


6.04 


0.0 


Ar +0.5 (s 2 p 5.5) 


13.82 


13.93 


0.8 


9.01 


9.02 


0.2 


Si+ (aV) 


31.61 


31.82 


0.7 


19.11 


19.19 


0.4 


Ga +0.5 ( s 2.0 p 0.5) 


46.97 


45.12 


3.9 


29.60 


29.97 


1.2 


Ge+ (s 2 ^ 1 ) 


30.38 


29.60 


2.6 


18.93 


19.12 


0.9 


Ti+ (sVrf 1 ) 


90.50 


90.61 


0.1 


47.98 


48.26 


0.6 


Quadrupole 














K+ 


18.7 


— 


— 


18.2 


— 


— 


K +0.3 (s 0.7) 


2530 


2521 


0.3 


3044 


3029 


0.5 


C+ (sV) 


16.4 


16.4 


0.0 


16.8 


16.8 


0.0 


Ar +0.5 (s 2 p 5.5) 


38.2 


38.2 


0.0 


37.8 


37.8 


0.0 


Si+ (sV) 


107.5 


107.5 


0.0 


108.2 


108.3 


0.1 


Ga +0.5 (s 2.0 p 0. 5) 


220.1 


218.8 


0.6 


228.5 


228.1 


0.2 


Ge+ (sV) 


107.4 


106.7 


0.7 


108.9 


108.5 


0.3 


Ti+ (sVrf 1 ) 


301.8 


286.8 


5.0 


319.9 


316.1 


1.2 
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FIGURES 

FIG. 1. Magnitudes of some important hardness matrix elements for Ar in configuration 
3s 12 3p 57 3d 01 . The hollow bars represent the all-electron results while the dashed bars repre- 
sent pseudopotential results at r cs = 0.9 a.u., r cp = r cc i = 1.0 a.u. Different contributions are 
decomposed into different columns as, "h" (Hartree), "xc" (exchange-correlation), "FWA" (total 
hardness under frozen wavefunction approximation), "SC" (self-consistent correction), and "total" 
hardness. The L > hardness matrix elements are rescaled by a factor of 2L + 1 to make them 
more clear. 

FIG. 2. Hardness matrix elements for Si + in configuration s 2 p 1 d°. Hollow bars represent AE 
results, while the dashed and shaded bars represent PSP results at r cs = 1.1 a.u., r cp = 1.2 a.u., 
fed = 0.8 a.u., and r cs = 2.3 a.u., r cp = 2.4 a.u., r c d = 0.8 a.u., respectively. 

FIG. 3. Hardness matrix elements for O in configuration s 1A p 35 d 01 . Hollow bars represent AE 
results, while the dashed and shaded bars represent PSP results at r cs = 0.8 a.u., r C p — 0.9 a.u., 
fed = 0.8 a.u., and r cs = 2.9 a.u., r cp = 3.0 a.u., r c d = 0.8 a.u., respectively. 

FIG. 4. Hardness matrix elements for K in configuration s - 25 ^- 25 ^ - 25 . Hollow bars represent 
AE results, while the dashed and shaded bars represent PSP results (r cs = 1.8 a.u., r cp = 2.3 a.u., 
r C( [ = 1.2 a.u.) without and with the LFC partial-core correction, respectively. 

FIG. 5. Hardness matrix elements for Ti in configuration s 1 p°' 25 d 2 . Hollow bars represent AE 
results, while the dashed and shaded bars represent PSP results (r cs = 1.8 a.u., r cp = 2.3 a.u., 
fed = 0.8 a.u.) without and with the LFC partial-core correction, respectively. 

FIG. 6. Hardness matrix elements for C in configuration 2s 1 " 9 , 2p , 3d " 3 . Hollow bars repre- 
sent AE results, while the dashed and shaded bars represent PSP results (r cs = r cp = 0.6 a.u., 
fed = 0.8 a.u.) without and with the LFC partial-core correction, respectively. 
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FIG. 7. Hardness matrix elements for Ga in configuration s 2 p°- 5 d°. Hollow bars represent AE 
results, while the dashed and shaded bars represent PSP results (r cs = 0.8 a.u., v cs — 0.9 a.u., 
r cc i = 1.2 a.u.) without and with the LFC partial-core correction, respectively. 

FIG. 8. Calculated average hardness error for Ar, Si, Ga, and K as a function of core radius 
used in the PSP generation. The configurations used to calculate the average hardness errors are 
Ar (Ss 1 -^- 95 ^- 05 ), Si (Ss^p^Sd - 1 ), Ga (4s L2 V- 7 3d - 1 ), and K (4s a85 4p ai 3d°- 05 ). 

FIG. 9. Calculated average hardness error for C and O as a function of core radius used 
in the PSP generation. The configurations used to calculate the average hardness errors are C 
(2s L9 2p 2 -°3d - 1 ) and O (2s L5 2p 4 -°3d - 5 ). 
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Fig 1, "Chemical Hardness ..." by Fillippetti et al. 
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Fig 2, "Chemical Hardness ..." by Fillippetti et al. 
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Fig 3, "Chemical Hardness ..." by Fillippetti et al. 
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Fig 4, "Chemical Hardness ..." by Fillippetti et al. 
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Fig 5, "Chemical Hardness ..." by Fillippetti et al. 



h xc FWA SC total 



ss,ss 
L=0 



PP,PP 
L=0 



I 



I 



ii 



I 



I 



i 



0.6 
0.3 
0.0 
-0.3 



dd,dd 
L=0 



sp,sp 
L=1 

(x3) 

sp,pd 
L=1 

(x3) 

P^P 

(x5) 

dd,dd 
L=2 

(x5) 



I 



UZLX 



l—hTTl— I 



M I U4J 



TZKJ 



Fig 6, "Chemical Hardness ..." by Fillippetti et al. 
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Fig 7, "Chemical Hardness ..." by Fillippetti et al. 
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